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The quark-meson coupling model for a nucleon is considered. The model describes 
a nucleon as an MIT bag, in which quarks are coupled to scalar and vector mesons. A 
set of coupled equations for the quark and the meson fields are obtained and are solved 



OO . in a self-consistent way. It is shown that the mass of a nucleon as a dressed MIT bag 



interacting with a- and cj-meson fields significantly differs from the mass of a free MIT 



O ! bag. A few sets of model parameters are obtained so that the mass of a dressed MIT bag 

^ I becomes the nucleon mass. The results of our calculations imply that the self-energy of 

the bag in the quark-meson coupling model is significant and needs to be considered in 

doing the nuclear matter calculations. 

^ ■ PACS number: 24.85. +p, 21.65.+f, 12.39.Ba 



1 Introduction 

Quantum hadrodynamics |[T|] proved to be a very successful model in describing the bulk 
properties of nuclear matter as well as the properties of finite nuclei. The model is rather 
simple with only a few parameters, and yet it has been successfully applied to a vast 
number of nuclear matter and nuclear structure problems. On the other hand, this model 
has a fundamental shortcoming that the nucleon with a composite structure is treated 
as a point particle. Several years ago a model to remedy this problem was proposed 
by Guichon p. He proposed a quark-meson coupling (QMC) model and investigated 
the direct quark degrees of freedom in producing the nuclear saturation mechanism in 
nuclear matter. The model describes nuclear matter as non-overlapping MIT bags P, Q 
interacting through the self-consistent exchange of scalar and vector mesons in the mean- 
field approximation. In other words, the nucleon in nuclear matter is assumed to be a 
static spherical MIT bag in which quarks interact with the scalar (a)- and the vector 
(u;)-meson fields. The model was refined later to include the nucleon Fermi motion and 
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the center-of-inass corrections to the bag energy and was apphed to a number of 
problems [6 — 11]. 

In the MIT bag model p, Q the bag constant B and a phenomenological parameter 
Z are fixed such that the nucleon mass of 939 MeV is reproduced for some bag radius 
R. The bag constant B produces the pressure to make a bubble in the QCD vacuum. 
Z is to account for various corrections including the zero-point motion. In the QMC 
model the quark-cr coupling constant (g^) and the quark-cj coupling constant (g^) are 
introduced. Usually, B and Z are fixed first so that the mass of the free MIT bag becomes 
equivalent to the the nucleon mass of 939 MeV for a certain bag radius. B^^^'^ and Z-l'^'^'^ 
thus obtained are used in the QMC model calculations for nuclear matter problems, g^ 
and g^ are then determined so as to reproduce the binding energy per nucleon (B.E. = 
— 16 MeV) at the saturation density {p% = 0.17 fm~^). 

However, due to the interaction between the quarks and the a- and the tu-mesons it is 
expected that the mass of a single dressed MIT bag in free space, when calculated with 
B-^^'^^ and Z^^^^ described above, would no longer be the nucleon mass of 939 MeV. Such a 
deviation of the mass of the dressed bag in free space has been neglected in previous QMC 
calculations in choosing the model parameters. If the deviation in mass is significant, it is 
necessary to modify the parameters to get the correct nucleon mass before implementing 
them in the nuclear matter calculations. It is this change in the mass of a dressed MIT 
bag that we address ourselves to in this paper. In Section 2 we describe the QMC model 
for a nucleon. A brief review of the QMC model for nuclear matter is given in Section 3. 
Section 4 contains the results of the QMC model calculations for a nucleon and nuclear 
matter. The paper is summarized in Section 5. 

2 Quark-meson coupling model for a nucleon 

The Lagrangian density for the MIT bag in which quark fields ip are coupled to the a- 
and uj- fields may be written as 
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i ^^Yd,^ - {d,^)Y^) - ^ {gll.oj'' + (m, - gla)) ^-B 



where F^^, = d^oj^ — dyOJ^, 6^ is the step function for confining the quarks inside the bag, 
and As is the 5-function at the bag surface. We will neglect the isospin breaking and 
take irtq = {niu + vn^jl hereafter. For numerical calculations m^ will be taken to be 
zero. The mass of a (cu) is taken as 550 (783) MeV. From this Lagrangian density the 
equations of motion for quark field ■?/', sigma field a, and omega field lo^ follow. For t/^ we 
have 

Yl^{id^ - glu,) - K - 9»] ^0. = ^(1 - ^rn,)^As. (2) 

The left hand side of this equation gives us the equation for the quarks inside the bag, 
and the right hand side provides us with the linear boundary condition at the bag surface 



{r = R). The equations for a and u^ are, respectively, 

d^d^a + mla = gl^^e, (3) 

and 

d''F,^ + mlu;^ = gli,^^^e,. (4) 

If we consider only the ground state quarks in the static spherical MIT bag and keep 
only the time-component of cj^, Eq. (Q) becomes the following coupled linear differential 
equations for g{r) and /(r) inside the bag; 



dfjr) 

dr 
dg{r) 



2l^ + E-{r)g{r) (5) 

r J 

E+(r)/(r), (6) 



dr 

where g{r) and /(r) are the radial parts of the upper and the lower components of i{j, 
i.e., 

^(t,r) = e--^*/^f .i^W ^^, (7) 

where a is the Pauli spin matrix and Xq is the quark spinor. Also, 

E'-ir) = |-^>o(r) + K-^Mr)) 

E-{r) = |-^>o(r)-K-^Mr)) (8) 

with cr(r) and u;o(r) being the a- and the time component of the a;-fields. Since we are 
dealing with the ground state, Eqs. (^ are the equations in the radial coordinate r only. 
The linear boundary condition from Eq. (|^), when rewritten by using Eq. (|^), reads 

at the bag surface, r = R, and determines the eigenvalue Xg of the quarks. e° is then 
given by Jxf+{Rrnqf = Xq for m-g = 0. 

In the static spherical approximation Eqs. (H) and (|) are reduced to 

(V2-mXr) = -(7^(3p,Mi?-r) (10) 

and 

(V^ - ml)u,{T) = -gl{'ipB)e{R - r), (11) 

respectively, where ps{= V"^) and Pb{,= '4'loi^) are scalar and baryon densities, respec- 
tively, and 3 is multiplied by them to account for the sum over 3 quarks. Equations (|^), 
(I), (|T0|), and (0) constitute a set of coupled equations for ip, a, and uo, which need to 
be solved self-consistently. 



By solving these equations we can obtain the eigenvalue of the quarks and the energy 
(-Eiv) of the nucleon bag. Ej^ can be computed by using 



Em = J dh Too (12) 

with 
Too = (I^tV + i?) e{R -r)-\ (iVuj.f + {m^uj.f) + 1 [{Vaf + {m^af) , (13) 

where Sq is given by 

^ = 3^ = 3^ (14) 

x? R R R R 

with the sum over 3 quarks taken into account. Correcting for spurious center-of-mass 

motion in the bag, the mass of the nucleon bag at rest is taken to be & O 



M^ = sjE% - (p2_), (15) 

where (pirn.) = J2k=i{Pk) = 3(xg/-R)^. By minimizing M^ with respect to the bag radius 
R, we can get the nucleon mass and the bag radius. 



3 Quark-meson coupling model for nuclear matter 

Now, let us consider the QMC model for nuclear matter. Detailed formulations and 
justification of the model can be found in Refs. P, |^, |, 0, and here we only briefly 
sketch the model for later discussions. In addition to the static spherical approximation 
assumed in Section 2 we now further use the mean-fleld approximation. From Eq. (0) 
the quark fleld i/j inside the bag in the mean-fleld approximation may be written as 

t^d^ - g^oo-f' - (m, - g^a)] tp = 0, (16) 



where a and u are the mean-flelds for a- and cu-mesons, respectively. Eq. (|16D has the 
following simple solution. 



.(M)^A^e— (,^/f^/f^^,,)-^, (., 



where e^ = i7g + g^ujR 



^' ^% Af-' = 2R'f,ix)[n,in,-l) + Rm;/2]/x' 



\ng + Rm* 






with Vtq = Jx"^ + {R m*)2 and m* = niq — g^a. The boundary condition 

Mxr/R)=PgMxr/R) (19) 



at r = R gives us the x value. Then the energy of the bag can be expressed as 



and the effective mass of the bag in nuclear matter is taken as 
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Using the minimum condition 
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(20) 

(21) 
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we can determine the effective mass. 

The total energy per nucleon at nuclear matter density, p^-, including the Fermi 
motion of the nucleons, can be written as [|, ^ 



E, 
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tot 



d'kjM^^g + k2 + 



(3^, 



q\2 



-PN + 



m 



'^a^ 



where the spin-isospin degeneracy 7 is 4 for symmetric nuclear matter. 
The mean field uj created by the nucleons is given by |[l|, 0, |^, |^ 



u 



mt 



and the scalar mean field a is determined by the thermodynamic condition 



da 



which yields the self-consistency equation 

m^ [ZTrj-^ 



0, 



R, pn 






Mil + k2 
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Equation (EBf) is also to be solved consistently with Eq. (O 
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4 Results and Discussions 

We first present the results for a single nucleon. In presenting the results we take the 
cases in which the bag radius is 1 fm. In Fig. 1 the nucleon mass Mj:q is plotted as a 
function of the bag radius R. The solid curve represents the free MIT bag mass {Mjj^^) 
without a coupling of the quarks with the mesons; da = 9^ = 0- ^^ the calculations B^^^ 
and Z are taken to be 136.3 MeV and 1.153 0, respectively, which give us the minimum 
of MjJ^^ as 939 MeV at 1.0 fm. If we include the interaction of quarks with meson fields 
using the couphng constants gl = 5.605 and g^ = 1.152 chosen by Saito and Thomas[|], 
for instance, a solution to the field equations (||), (||), (p!Ol), and ([TT| ) does not exist. In 
Eq. (^ both dg{r)/dr and /(r) are negative (or zero) for all radii, and thus E~^{r) needs 
to be positive. However, due to the large attraction coming from the large g^ value in 
relative to the small g^ the eigenvalue e^ in Eq. (^ becomes relatively small, which 
results in E~^{r) being negative for this set of parameters; hence the non-existence of the 
solution. Therefore we show here results for other sets of parameters. 

First, we include only the coupling of the quarks with the a-meson; g^ = 0. The 
bag mass decreases due to the attractive cr-coupling as shown in Fig. 1 by the thin 
dashed and the thin dash-dotted curves, respectively, for g^ = 2 and 4. M^- is reduced 
to 930 and 903 MeV, and the bag radius also slightly decreases to 0.997 and 0.979 fm, 
respectively, for g^ = 2 and 4. On the other hand, when the quarks are allowed to couple 
only with the cj-field, i.e., gl = 0, M^ increases as plotted by the thick dashed and the 
thick dash-dotted curves, respectively, for g^ = 2 and 4. M^ is pushed up to 954 and 
994 MeV due to the repulsive coupling, and the bag radius becomes 1.00 and 1.03 fm, 
respectively. 

If the quarks are allowed to couple with both the a- and the a;-mesons, cancellation 
between the attraction and the repulsion takes place. In Fig. 2 we show M^v calculated 
with g1 = g^ = 2 by the dashed curve in comparison with MJJ^^ {gl = g^ = 0) plotted 
by the solid curve. Mn for g1 = g^ = 2 has a minimum at 946 MeV and R = 1.01 fm. 
The mass difference (946 — 939 = 7 MeV) is just about what is left after the cancellation 
between the thin dashed curve {g^ = 2, g^ = 0) and the thick dashed curve {g^ = 0, 
g^ = 2) in Fig. 1. Mjy calculated with g^ = g^ = 4 is also plotted in Fig. 2 by the 
dash-dotted curve. The minimum of M^v is 964 MeV at i? = 1.02 fm. Again, the mass 
shift (964 — 939 = 15 MeV) is just about the repulsion left after the cancellation between 
the two dash-dotted curves in Fig. 1. It is worth noting here that when we use the same 
value for g^ and g^ the net result is a repulsion. This may be understood by looking at 
the meson fields. In Fig. 3 cr(r) and ujQ{r) are plotted as a function of radius for a fixed 
bag radius i? = 1 f m when g^ = 9^^ = 2. The figure shows that the shapes of cr(r) and 
c<Jo(r) are largely dictated by those of the source densities ps and pb, respectively. Just 
as Pb is bigger than ps at larger radii, uJo{r) is greater than cr(r) at larger radii, which 
is the significant radial region. Thus, when the two coupling constants are equal, a net 
effect turns out to be repulsive. 

In the nuclear matter at normal nuclear densities the internucleon distance is about 
2 fm. In the mean-field approximation meson fields are assumed to be constant in the 
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matter. On the other hand, Fig. 3 (a) shows that the range of the meson fields is 
roughly around 1.5 fm for both mesons. The fields become almost negligible at r ^ 2 
fm. However, since the nucleons in nuclear matter are in the Fermi motion, the meson 
fields with the range of about 1.5 fm may be approximated to mean-fields in the matter. 
Therefore, the present study of a single nucleon in the QMC model seems compatible 
with the QMC model for the nuclear matter, which will be discussed shortly. 

We then looked for parameters that could give us the correct nucleon mass for R = 1 
fm. First we fixed B^^'^ and Z as 136.3 MeV and 1.153 0, respectively, and searched 
for g'^ and g^. It was not possible to find g'^ and g^ that would yield the correct nucleon 
mass with the given values of B^^^ and Z. When we let Z also vary, keeping B^^^ fixed 
as 136.3 MeV, it was still impossible to find parameters producing the correct nucleon 
mass. On the other hand, when we let B^^^ vary, keeping Z as 1.153, we could obtain 
the correct mass with the parameter set 2 listed in Table 1. There are four parameters 
to be fixed, whereas there are only two quantities to fit; the nucleon mass and the bag 
radius. Thus the parameters are not unique. We list in Table 1 other sets of parameters 
that can produce the nucleon mass of 939 MeV for R = 1 fm. Mat's calculated with the 
parameter sets 2-7 are all indistinguishable from the solid curve in Fig. 1 or Fig. 2, 
so they are not plotted here. One can see from Table 1 that g'^ is not much larger than 
g^. g1 is often only slightly bigger than g'^, as the relativistic Dirac phenomenolgy has 
shown. In fact, this is expected from Fig. 2, which shows that with the equal values of g^ 
and g'^ there is only a small repulsion. To restore the correct mass, g^. would have to be 
a little larger than g'^. This is in contrast to the values of g^ and g'^ previously obtained 
by QMC model calculations for nuclear matter. Too large gl values compared to small 
glj were one of the concerns in the previous results of the QMC model calculations. 

Similar studies were made at other bag radii ranging from 0.6 to 1.0 fm, and similar 
results were obtained. Our results show that the change in the nucleon mass due to 
the coupling of the quarks with the mesons is not negligible. Such self-energy effects on 
the nucleon mass need to be taken into consideration in choosing the parameters for the 
calculation of nuclear matter properties. Thus, we performed the nuclear matter calcu- 
lations as sketched in Section 3 along with the calculation for a nucleon. As mentioned 
earlier, when we use g^ and g^ taken from Ref. [§, a solution for a nucleon problem does 
not exist. On the other hand, if we use the parameter sets 2 - 7 in Table 1, the calculated 
binding energies are not correct. Even saturation does not occur for these parameter 
sets. Therefore, we attempted to look for a new set of parameters that could reproduce 
simultaneously the nucleon mass and the nuclear binding energy. We found that it was 
not possible to get such parameters, at least, within the present scheme. Therefore, this 
additional requirement of reproducing a single nucleon mass can be used as a guide in 
refining the QMC model for nuclear matter. 



5 Summary 

Recently the QMC model has been frequently used to describe the properties of nuclear 
matter and finite nuclei. However, in the previous calculations the change in the mass 
of the bag due to the self-energy has been ignored. Thus we have applied the QMC 
model to a single nucleon. Our calculations suggest that this change in the nucleon mass 
is not negligible. Therefore, we performed the QMC model calculations not only for a 
nucleon but also for nuclear matter and looked for parameters that would give us both 
the nucleon mass and the nuclear binding energy correctly. However, we only found that 
it was not possible to obtain such parameters within the present framework of the model. 
Taking into account this change in the nucleon mass may provide us with an additional 
constraint on the model. 
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Table 1: The parameter set 1 is taken from Ref. [6]. The paremeter sets 2-7 are 
the parameters that can produce M^ = 939 MeV at the bag radius R = 1 fm. These 
parameters are not unique, so here we hst some of them. 



set 


5^4 (MeV) 


Z 


9l 




Mn (MeV) 


1 


136.3 


1.153 


5.605 


1.152 


No solution. 


2 


136.5 


1.153 


3.396 


2.490 


939 


3 


136.5 


1.169 


3.300 


2.616 


939 


4 


136.5 


1.183 


3.199 


2.723 


939 


5 


136.5 


1.199 


2.941 


2.760 


939 


6 


136.5 


1.170 


1.965 


1.752 


939 


7 


136.4 


1.238 


3.840 


3.746 


939 
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Figure 1: The nucleoli mass (M^v) calculated by the QMC model for a nucleon is plotted 
as a function of the bag radius R for different coupling constants g^. and g'^. B^/'^ = 136.3 
MeV and Z = 1.153 are used. 



Figure 2: The nucleon mass {Mi\i) calculated by the QMC model for a nucleon is plotted 
as a function of the bag radius R for different coupling constants g^ and g^. B^^'^ = 136.3 
MeV and Z = 1.153 are used. 



Figure 3: (a) The a- and the cu- meson fields calculated by the QMC model for a nucleon 
are plotted as a function of radius r for bag radius i? = 1 fm with g'^ = g^ = 2, 
^1/4 _ ]^3g_3 jvIgV and Z = 1.153. (b) The scalar and the baryon densities are plotted 
by the dashed and the solid curves, respectively, with the same parameters. 
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